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This paper presents a technique for investigating the stability and damping 
present in centrally preloaded radially symmetric multi-mass flexible rotor bearing 
systems. In general, one needs to find the eigenvalues of the linearized pertur- 
bation equations, though zero frequency stability maps may be found by solving as 
many simultaneous non-linear equations as there are dampers; and in the case of a 
single damper, such maps may be found directly, regardless of the number of degrees 
of freedom. The technique is illustrated for a simple symmetric four degree of 
freedom flexible rotor with an unpressurized damper. This example shows that whereas 
zero frequency stability maps are likely to prove to be a simple way to delineate 
multiple solution possibilities, they do not provide full stability information. 
Further, particularly for low bearing parameters, the introduction of an unpressur- 
ized squeeze film damper may promote instability in an otherwise stable system. 


INTRODUCTION 

The use of centrally preloaded squeeze film dampers for the attenuation of the 
unbalance response in turbomachinery has been well documented, and solution tech- 
niques which enable all equilibrium operation possibilities to be conveniently por- 
trayed for general multi-degree of freedom rotor bearing systems are increasingly 
available (refs. 1, 2). However, the question as to which of these equilibrium 
solutions is stable has not been as fully addressed. Indeed, earlier stability 
investigations for simpler squeeze film damped flexible rotors (ref. 3) showed that 
instability (in the linear sense) was indeed possible with unpressurized dampers 
below and above the first pin pin critical speed, though no instability was noted for 
pressurized dampers with retainer springs. The utility of squeeze film dampers to 
accommodate the influence of gyroscopic effects, non-rigid bearing mounts and super- 
critical operation on stability needs to be better quantified. Hence, it is the 
purpose of this paper to present a straightforward but general technique for investi- 
gating the stability and degree of damping present in general multi-mass flexible 
rotor bearing systems incorporating one or more centrally preloaded squeeze film 
dampers. The technique will be illustrated for a simple symmetric four degree of 
freedom flexible rotor. 



SYMBOLS 


A,B 

square matrices defined by 

A 

c 

= T -1 CT* 


equations (27) and (28) 

c* 

matrix defined by equation (18) 

C 

radial clearance of damper 

Cl 

= Cj/CCnij +m 2 )aj] 

c 

matrix of viscous damping and 
gyroscopic coefficients 

2c l 

viscous damping of disc in 
figure 3 
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etc 


c rs> 

etc. 

“ Qj^gi^/E(m^ "f"m 2 )(A3]j etc* 

e i 


eccentricity of i^ damper; 

1 } • • * j m 

*1,2 


unbalance excitation forces 



defined by equation (31) 


Ir 

Vectors of system excitation 
and hydrodynamic forces in 
stationary and rotating frames 
respectively 


6 Ir 

vectors of perturbated system 
excitation and hydrodynamic 
forces in stationary and 
rotating frames respectively 

F Si’ 


elements of Fg and F R 
respectively; i=l , . . . ,n 

i 


subscript (omitted where 
meaning is clear) 

j 


subscript or depending on 

context 

K 


system stiffness matrix 

K* 


matrix defined by equation (19) 

K l, 2 


* k x 2 /C (m 1 +m 2 )o> 2 ] 

K rs> 

etc. 

= P R34/t +m 2^ a)2 ^ etC * 



rotor stiffness in figure 3 

k 2 


retainer spring stiffness in 
figure 3 

L 


damper length 

M 


system mass matrix 

A 

M 


= T _1 MT* 

M 

1 -1,2 


= +m 2 ] 

m 


number of hydrodynamic damper 
stations or some character- 
is tic system mass 

m 


matrix defined by equation 
(All) 

A 

m 


matrix defined by equation 
(A12) 

2m x 


lumped mass of disc in figure 3 

m 2 


lumped mass at bearing stations 
in figure 3 

n 


number of system degrees of 
freedom (necessarily even) 


P S» P R 

matrix of damper stiffness 
coefficients in stationary and 
rotating frames respectively 

P Sij’ p Rij 

elements of Pg and P R defined 
by equations (7) and (20) 
respectively 

<*s> Qr 

matrix of damper damping coef- 
ficients in stationary and 
rotating frames respectively 

Qs ij’ QR ij 

elements of Qg and Q R defined 
by equations (8) and (21) 
respectively 

R 

radius of damper journal 

S, R 

subscripts denoting stationary 
(XYZ) and rotating (xyZ) 
frames respectively 

r, s 

subscripts denoting damper 
degrees of freedom x 3 and x 4 
respectively 

T 

transformation matrix defined 
by equation (11) 

x* 

matrix defined by equations 
(A4) to (A6) * 

t 

real time 

u 

unbalance parameter 
= p 1 m 1 /C(m 1 +m 2 )C] 

u 

state vector defined by 
equation (26) 

XYZ, xyZ 

stationary and rotating 
cartesian reference frames 
respectively 

X, x 

system displacement vector in 
stationary and rotating frames 
respectively 

X . x 
-o* -o 

system steady state displace- 
ment vector in stationary and 
rotating frames respectively 

6X, 6x 

perturbed system displacement 
vector in the stationary and 
rotating frames respectively 

x i 

elements of X or x 

x x , etc. 

= x x /C, etc. 


phase difference between i tn 
lumped mass and the stationary 


frame as defined in figure 2 
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£ damper eccentricity ratio = e/C 

C damping ratio at rotor mid-span 

in figure 3 - c 1 /(2m 1 to c ) 

Q eigenvector of the perturbed 

system as defined by equation 
(29) 

A transformation matrix as 

defined by equation (12) 

A* matrix defined by equation (A5) 

X system eigenvalue 

y lubricant viscosity 

p 1 disc mass eccentricity in 

figure 3 

T non-dimensional time = 0)t 

<j> angle between axes OX and Ox 

as defined in figure 2 

X eigenvector of the perturbed 


system as defined by equation 

( 22 ) 

$ Phase difference between 

rotating and stationary frames 
as defined in figure 2 

Q natural frequency at stability 

threshold 

0) rotor speed 

U) c a characteri stic sy stem 

frequency = /k 1 /m 1 in figure 3 

a bearing parameter 
= yRL 3 /C(m 1 +m 2 )C 3 ] in figure 3 

0) r - /k 2 / (m x + m 2 ) 

• denotes differentiation with 

respect to time t 

1 denotes differentiation with 

respect to time T 


THEORY 

Figure 1 depicts a general n degree of freedom rotor bearing system incorpor- 
ating one or more squeeze film dampers. For the system, the equations of motion, in 
fixed cartesian co-ordinates, may be written as: 

MX + CX 4* KX - F s (X,X) . (1) 

Assuming (i) axially symmetric rotor and foundation stiffnesses, (ii) viscous 
damping, (iii) central preloading of the hydrodynamic dampers, (iv) synchronous 
unbalance excitation, (v) negligible torsional and axial vibration, one can obtain 
circular synchronous solutions Xq to equation (1) as explained in reference 2, i.e. 

MXq + CXo 4* KX 0 = I S (?o»io> ^ ( 2 ) 


Note that X Q is a function of time. It is the degree of stability of these steady 
state solutions Xq that is of interest. 


If the steady state solution X Q is perturbed by 6X to X, whereupon Fg(X 0 ,X 0 ) 
changes to F g (X,X), i.e. if: 


X = X c 4- 6X , (3) 

then substitution of eqtxation (3) in equation (1) and utilizing equation (2) yields: 

MSX + CSX 4* KSX - 6f s , (4) 

where 6F S - F S (X,X) - 2 S (X 0 ,X 0 > (5) 

= P s SX + Qg6X + higher order terms , (6) 
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Figure 1 Typical multimass flexible rotor running in damped flexible supports 

3F S . 

and where P s _ = , (7) 

9F S* 

and Qc . . » — ^ . (8) 

Thus, neglecting higher order terms, equation (4) may also be written as: 

MfiX + (C -Qg)6X + (K-P s )6X = 0 . (9) 

The existence of the partial derivatives is assumed and they are evaluated at X 0 ,X Q . 

Unfortunately, the elements of Qg and Pg, involving derivatives of damper 
forces, are in general time dependent, so that equations (9) do not reduce to an 
eigenvalue problem. To overcome this difficulty, one may choose a rotating cartesian 
reference frame (x,y,Z) wherein the cartesian axis pair (x,y) rotates with, angular 
velocity U) about the Z axis as shown in figure 1. If x be the vector of displace- 
ments in the rotating frame, then: 

X = Tx , (10) 
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If x Q be the steady state solutions of equation (15), then following the same 
arguments as above for the stationary frame, the equations of motion consequent upon 
system perturbation 6x, become: 

M6x + (C + 2wM - Q r )6x + (-w 2 M + ojG + K - P R ) fix = 0 , (16) 


or 

where 


and 


MSx + C*6x 4 - K*6x = 0 , 

C* = C 4 - 2a)M - Q r , 
K* = K + coC - co 2 M - P R , 




9F- 


R-i 


3x^ 


y 




(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


Again, the P R . . and the Qr. . are assumed to exist and are evaluated at x . (Note 
that x Q = 0 . ) 


In general, the steady state damper eccentricities e^ will have constant phase 
angles (y^ - 1 p) with respect to the x axis, as shown in figure 2. Hence, the elements 
of P R and Qr will be functions of e± and (y^ - ip) , where \p is arbitrary. By selecting 
\p equal to one of the Yj/s (say Yi)> Pr and Q R will, in the general case of m 
dampers, be functions of (2m - 1) quantities e 2 ,e 2 , . . . ,e m , Y 2 “Y 1 * >YnrYi • 

Note that the coefficients of 6x, 6x and 5x in equation (17) are constants, 
i.e. independent of time, so one is immediately in a position to investigate the 
stability and the damping pertaining to the equilibrium solutions x Q by examining the 
solutions of equation (17), a set of n homogeneous second order linear differential 
equations with constant coefficients. Thus, by assuming solutions of the form: 

Sx = , (22) 

where the X and A may be complex, non-trivial solutions of equation (17) exist, if, 
and only if: 

detCA 2 H + AC* + K*] = 0 . (23) 

Equation (23) , the characteristic equation of the perturbed system, is a poly- 
nomial of degree 2n in A. The stability of, and the damping pertaining to the 
equilibrium solutions x Q depend on the real parts of the roots of equation (23) with 
the stability threshold being determined by that combination of system parameters 
which result in any pair of roots, A x 2 = ±jf2, where 9, is the natural frequency at the 
stability threshold. Note that one such possible stability threshold corresponding 
to 9-0 (i.e. A x 2 = ±j0) will occur whenever: 

det[K*] = 0 . (24) 


If stability threshold determinations were the sole requirement, one could 
dispense with the need to find the roots A 1 ,...,A 2n and apply a linear systems theory 
technique, e.g. Routh’s Criterion, to the coefficients of equation (23). Such was 


396 



the procedure adopted in reference 3. If the degree of damping is of interest, one 
needs to resort in general to some numerical procedure for finding the roots. 

An alternative and numerically far more simply approach, in so far as program- 
ming effort is concerned, is to utilize generally available eigenvalue solution 
procedures, by recasting equations (17) into a recognizable eigenvalue problem. The 
normal procedure is to transform these n second order differential equations into the 
2n first order equations (see for example reference 4) : 


Au + Bu - 0 


(25) 


where 


u 


<5x 

6x 


(26) 


A 


0 M 
M C* 


(27) 


and 


B = 


By assuming solutions of the form: 

u = 


-M 0 

. 0 K* 



(28) 


(29) 


where the i] and X may be complex, non-trivial solutions of equation (23) exist if, 
and only if: 


det[XA + B] = 0 . (30) 

Equation (30) is a polynomial of degree 2n in X, It is equivalent to the 
characteristic equation (23). The 2n values of X, the roots of the characteristic 
equation, are in this formulation more commonly referred to as eigenvalues. Once the 
2n eigenvalues have been located for a given choice of system parameters, both the 
stability and the degree of damping present can be quantified. 


SYSTEM WITH ONE DAMPER 

So far the problem formulation has been quite general. If one restricts 
attention to a single damper, an important simplification occurs if ip is set equal to 
Yj, whereupon the elements of P R and Q R in equations (16) are functions of e x only. 
Thus, the need to evaluate is avoided, thereby simplifying stability evaluation 
and design data portrayal. 

For design study purposes, it is convenient to non-dimensionalize equations 
(16) by dividing the force equations by mCW 2 and the moment equations by mC 2 w 2 and 
introducing the non-dimensional time x=CJt. The above theory is otherwise unchanged, 
except all quantities are now non-dimensional. In particular, one can utilize 
equation (24) immediately to find those stability thresholds which pertain to zero 


397 



natural frequency. Note that the non-zero elements of are located in a 2x2 sub- 
matrix, the elements of which are products of 0)^/0) and functions of e, where 0)^ is a 
bearing parameter involving the bearing dimensions and lubricant viscosity. Thus, 
whereas equation (24) is a non-linear equation in £, and for a given 0)^, could be 
solved iteratively by some appropriate technique such as the search procedure, one 
may note that for some assumed e but unspecified (%, this equation is a quartic in tofo 
when the damper is flexibly supported, and a quadratic in when the damper is 
rigidly supported. Hence one can determine directly from equation (24) the physical- 
ly meaningful values of oo^ (if any) which satisfy it. A repetition of this for other 
values of U) would enable the stability thresholds corresponding to zero natural 
frequency and various values of E to be drawn on a design map with 0)^ as ordinate and 
(0 as abscissa. 

Note that such a map is possible regardless of the number of degrees of free- 
dom, or of the unbalance distribution. If drawn to the same scale, it can be direct- 
ly superimposed on the corresponding equilibrium orbit eccentricity design maps in 
reference 2, thereby indicating at a glance, the likelihood of operation in the 
vicinity of, or within, this particular type of unstable region. Note that the 
absence of other stability thresholds (at some non-zero natural frequency) has not 
been proven, and to be sure that all stability regions have been located, one would 
need to apply a more general technique, e.g. Routh's Criterion, to equation (23) over 
the range of 0 < e < I, for the range of values of and 0) of interest, a rather 
daunting task even for systems with only four degrees of freedom. 


ILLUSTRATIVE EXAMPLE 

The utility of the above approach will be illustrated for the single disc 
symmetric flexible rotor, previously investigated for stability in reference 3 and 
for which equilibrium solutions are available in reference 2. 

Figure 3 is a diagram of this system with node 1 being taken at the central 
disc of mass 2m 1 , and nodes 2 and 3 at the ends of the rotor which are supported by 
identical squeeze film dampers. The lumped mass at the bearing ends is m 2 , the 
retainer spring for central preloading has stiffness k 2 and the rotor stiffness 
between the central and either end node is k a . All unbalance is assumed to be at the 
disc, resulting in a disc mass eccentricity p x . Viscous damping at the disc is 2c j. 
Since the rotor is symmetric about the disc, it will suffice to consider one half of 
the System only. The equations of motion, appropriately ordered, for this systein at 
a rotor speed U) are given by: 

m l X 1 + Cl X x + k 1 (X 1 - X 2 ) 

*4*2 + c x X 2 + k x (X 2 - x 4 ) 

m 2 X 3 + k x (X 3 - X x ) -f k 2 X 3 

m 2 X 4 + k x (X^ - X 2 ) + k 2 X 4 

The equations of motion for the perturbed system are given by equations (16), 
which, in nondimens ional form, upon dividing each equation by (m x + m 2 ) CU) 2 and letting 
T — cot , become: 


p 1 m 1 0) 2 cos(a)t) 

p 1 m 1 a3 2 sin(a)t) 

F. 




= F 0 


(31) 


= F, 
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Figure 3 Flexible, symmetric unbalanced rotor supported on identical 
squeeze film dampers and retainer springs 



where, to simplify the notation, the damper degrees of freedom x 3 and are denoted 
by r and s respectively* The coefficients K rr , C rr , . etc. are available in 
analytical form for unpressurized and fully pressurized (tt and 2tt film) dampers for 
both the short (Ocvirk) and finite width (Warner) bearing approximations (ref. ,3). 
Thus, for the unpressurized or fully pressurized bearings respectively, and using the 
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short bearing approximation, one has: 


„ 4e(l + e 2 ) 

^rr ~ (1 - e 2 ) 3 ° r 

0 , 

(33) 


TT 

(34) 

K rs “ U SS “ 2(1 - £2)3/2 or 

(1 - £2)3/2 > 

„ „ — TT ( 1 + 2 e 2 ) 

^sr c rr 2(1 - e 2 ) 5 / 2 ° r 

— tt ( 1 - 2 e 2 ) 

(1 - e 2 )V 2 » 

(35) 

— 2e 

^ss *“ ^rs “ ^sr ' (i - e 2 ) 2 or 

0 . 

(36) 


The equation obtained by setting the determinant of K* in equation (32) equal 
to zero is a quadratic in (%, and may be solved for given values of £ and 0 ). 

To present the data in terms of non-dimensional quantities, it is convenient 
to plot a non-dimensional bearing parameter against the non-dimensional rotor 

speed U)/u> c , where U) c is some characteristic natural frequency of the system, say, the 
highest undamped natural frequency of the system which will fall below the desired 
operating speed. Equation (24) was then solved for o)|j/cl) g over a range of U)/(jO c and £ 
to obtain the zero natural frequency stability threshold maps if figure 4, using the 
following values of the non-dimensional system parameters: 

M 2 = 0.25 , 

M x = 1 - M 2 = 0.75 » 

Kj = (1 - M' 2 )/(o>/o ) c ) 2 = 0.75(to/w c ) 2 > 

K 2 = (u r /w c ) 2 /(w/w c ) 2 = 0 . 25 (co/co q ) 2 » 

Ci = 2^(1 - M 2 )/(w/w c ) = 0.0075/(w/w c ) • 

The quantities w r , w c and C are defined in the notation. (Using the 
notation in reference 3, the above choice of parameters would correspond to a =0.25, 
f -0.5, C = 0.0005, a= 0 )/a) c and B = 005 /co c . ) Here, co c is the first pin-pin critical 
speed of the rotor. 

Note that provided the same scales are used for the axes ooj > /oo c and Q)/u> c , one 
can overlay this stability threshold map not only over the equilibrium eccentricity 
orbit map given in figure 3 in reference 2 which is for an unbalance parameter U = 0.3, 
but also over any such map for the system, regardless of the unbalance distribution. 

* These zero frequency stability thresholds were compared to the stability thresholds 
for o)b/u> c = 0.1, 0.3, 0.6 and 1.0 in figure 8 of reference 3, which were determined 
using Routh’s Criterion. No stability thresholds were found for £<0.61. For U3 /(jl) c <1 
there is agreement. However, for U)/03 c > 1, the predicted regions of instability in 
figure 8 of reference 3 exceed those of figure 4. Thus, figure 4 does not provide 
full stability threshold information; and for this, one would need to resort to an 
approach such as Routh f s Criterion. 

As an alternative or in addition to seeking stability thresholds, and to 
determine the degree of damping pertaining to the equilibrium solutions, one can find 
the eigenvalues A-j^ for the system in figure 3 by forming equation (25) and using an 
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Figure 4 Zero natural frequency stability map for the system in figure 3 


eigenvalue solver to determine the X± (±-l 9 2, * . . , 8)* This was done for the 
particular equilibrium solutions indicated on figure 5 (corre spending to figure 4 in 
reference 2). If any of the X± have non-negative real parts, the solution is 
regarded as unstable. Where all the X± have negative real parts the solution is 
stable* Of particular interest is the degree of relative damping present in such 
stable solutions, particularly in the low orbit eccentricity solution when multiple 
solutions are possible, as this gives a qualitative indication of the likelihood of 
jumping to the undesirable and possibly unstable high orbit solution upon some system 
perturbation. Various means of determining this relative damping are available (ref .4) 
such as the smallest logarithmic decrement, -27rRe(A) / | Im(A) | , or the smallest damping 
ratio, -Re(A)/|A|. Either of these quantities requires the determination of the A^. 

The latter quantity xs indicated for illustrative purposes in figure 5. For example, 
at a)/oo c =0*3, the damping ratios corresponding to ca^/a) c =0.006, 0.06 and 0.3 are 
0.0053, 0*047 and 0.23 respectively. The equilibrium solution for (A>k/u> c =0.3 could 
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Figure 5 Damping ratios at selected points on the rotor response curves for 
various values of bearing parameter for the system in figure 3 


be regarded as having a high degree of relative damping whereas that for 0%/(j0 e “0.006 
could be regarded as having a low degree of relative damping. 

Wherever multiple solution possibilities exist, all intermediate solutions 
were unstable. In all such cases, the oscillation frequency was zero, and the 
instability was predicted by the stability map of figure 4. However, unstable 
single solutions are also indicated, as well as unstable higher eccentricity 
solutions in case of multiple solutions. In such cases the introduction of the 
damper has worsened system behaviour. Such unexpected unstable solutions occurred 
for a )b/u) c equal to 0.006 and 0.06 but not for 0.3. In all such cases, the 
oscillation frequency was always non-zero and was not predicted by the zero frequency 
stability map of figure 4. 
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This provides additional proof that figure 4 does not provide full stability 
threshold information. The number of equilibrium solutions investigated for 
instability proved insufficient to yield definitive trends in the variation of 
relative damping as one progresses along any particular frequency response curve in 
figure 5. To pick up such trends, a more thorough investigation of the would 
need to be undertaken. This is not warranted, for the purpose of this example was 
not to investigate in detail the particular system of figure 3, but rather to 
illustrate how simply the technique developed in the paper may be put to practical 
use. 


Note that although figure 4 does not represent global stability thresholds, 
it does predict the instability of all the "intermediate" solutions and so, it 
serves as a locator of jump speeds, i.e. speeds at which there is a transition from 
a speed for which there is only one solution to one where there are at most two 
stable solutions or vice versa. Since the high orbit eccentricity solution, if 
stable, is still undesirable, generally resulting in unbalance force magnification, 
operation in the vicinity of such jump speed regions should be avoided. Though only 
proven for the model in figure 3, this equivalence between the zero frequency 
stability map and jump speed location is expected to be valid for general multi- 
degree of freedom systems. If so, a relatively simple way has been found for 
delineating multistable operation possibilities for any system with one damper. 


CONCLUSIONS 

1. A technique is developed for investigating the stability of and the 
degree of damping in the circular synchronous orbit solutions of n 
degree of freedom rotor bearing systems. In general, the technique 
requires finding the 2n eigenvalues of the linearized perturbation 
equations. 

2. The perturbation equations of motion are not a function of the 
unbalance distribution, so for a given system a single global 
stability map suffices for all unbalance distributions of interest. 

3. Zero-frequency stability thresholds may be found by solving as many 
simultaneous non-linear equations as there are dampers. If the system 
contains one damper only, all such stability thresholds may be found 
directly by solving at most a quartic equation. 

4. Zero-frequency stability maps do not provide full stability 
information, but for the four degree of freedom system investigated 
in the illustrative example, and probably for more general higher 
degree of freedom systems as well, such maps provide a simple way to 
delineate multiple solution possibilities. 

5. Depending on the system parameters, single equilibrium and, where 
multiple equilibrium solution possibilities exist, the high orbit 
eccentricity solution may also be unstable, with the likelihood of 
instability apparently increasing as the bearing parameter is reduced. 
Thus, the introduction of an unpressurized squeeze film damper may 
promote instability in an otherwise stable system. 
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APPENDIX 


Since 


where 



T = 


Tx , 


Tx + Tx , 

Tx + 2Tx + Tx , 


A* v . 


V 


0) 


A*. 


9 




L o — 

— A* J 

and 

A* = 

-sin (j) 
cos (j) 

-coscj) 
-sincj) _ 

it follows that 

T = goT* , 


and 

T = -ga 2 T . 



Hence, substitution of equation (10) into equation (14) gives: 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 
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(A8) 


T -1 MTx + (2cJI -1 Mr* +T _1 CT)x + (uT _1 Cr* +T~ 1 KT -a) 2 T -1 MT)x = Fr , 


where 



(A9) 


and 


a- 1 


cos(f> sin^) 
-sin <J) coscj) 


(A10) 


To evaluate T^MT, it is convenient to partition it into 2x2 submatrices. 

Any one of these submatrices will have the form A^mA, where m is the corresponding 
submatrix of M. Then it is easy to show that A^mA equals m, if, and only if, the 
elements of m are of the form; 


mi -m2 
m2 mi 


(All) 


Hence, if equation (All) is satisfied for all submatrices of M, then T _1 MT 
equals M. Conditions similar to equation (All) are required for T^CT and T-^KT to 
equal C and K respectively. As may be seen from reference 6., such conditions are 
satisfied for M, C and K matrices in general, even when gyroscopic effects are 
present . 


To evaluate T~ 1 MT* , it is again convenient to partition it into 2x2 sub- 
matrices. Any one of these submatrices will have the form A^mA*, where m is the 
corresponding submatrix of M. Hence, noting that m will be of the form given by 
equation (All), A"~ l mA * will equal m where m is given by: 


A 

m 


—m2 -mi 

mi -m 2 


(A12) 


__ , /V A 

Thus T MT* equals M, where the 2x2 submatrices m of M are formed from the 
corresponding submatrices 1 of M according to equations (All) and (A12) . 

Similarly T“ T CT* equals C. Hence, equation (A8) simplifies to: 

Mx + (2aM +C)x + (-0) 2 M + 10 C +K)x = Fr . CA13) 

Equation (A13) is the equation of motion in the rotating frame. Note that the 
coefficients of x, x and x are all constants. 
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